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The design and analysis of CPV systems require knowledge of direct normal solar irradiation but 
ground-based measurements of these data are only available for very few locations. Nowadays, 
meteorological databases that estimate direct normal irradiation from satellite images and other data 
sources are used. However, values provided by the different existing databases show large dispersion 
due to different estimation methods, input data and base years. In this paper, a simplified method for 
calculating direct normal irradiation is presented. It has been obtained from previous models proposed 
by several authors. One of its advantages is that it only requires latitude and global horizontal 
irradiation as input data. As global irradiation is easy to find or measure, the procedure becomes a 
useful tool in renewable energy applications. The accuracy of this method is similar to that of the 
existing databases and it is able to easily generate a mass of direct normal irradiation data for different 


areas worldwide. 


© 2012 Elsevier Ltd. All rights reserved. 


Contents 


Qu WN 


Description of the proposed method..................0.000 00s 
Results and discussion. .... 0... 00.0 cece eee eee eee 
Application: eris tenenan ana EE E ben ee eds hada tage 
CONCIUSIONS: ss sepi cae eae dain e a Bede ee a DREGE 
Acknowledgments .......... 00. c cece cece cece eee n ene nees 
References: sss eker nels d dia a Hie EE nate 14 ete anaes 


1. Introduction 


Photovoltaic devices are at a mature technological stage with a 
long and wide field experience, but if this energy source wants to 
compete against other renewable energy sources or even to get 
into the traditional energy generation system, further scientific 
progress in the behavior of these kinds of systems is necessary. 
Concentrator Photovoltaic (CPV) technology can have a main role 


Abbreviations: CPV, concentrator photovoltaic; DNI, direct normal irradiance or 
irradiation; TMY, typical meteorological year; HC, HelioClim; HCPV, high con- 
centration photovoltaic 
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in this challenge, as it has proved, in recent published research, to 
have the potential to achieve high levels of energy conversion 
performance [1]. 

CPV is based on the use of optical devices that increase the 
light received on the solar cell surface. The concentration 
mechanism is achieved by lenses or mirrors that either refract 
or reflect the incoming light from the sun on top of the cell 
receiver surface. Depending on the concentration process, as well 
as on the concentration area, there is a wide range of optical 
device configurations that can be used in the implementation of a 
CPV module [2]. 

The use of concentrators implies that CPV systems only work 
with the Direct Normal Irradiance (DNI). So it is necessary to 
know DNI data in order to estimate the energy that will be 
produced by the system, perform economic analysis, supervise 
plant operation, etc. However, DNI ground-based measurements 
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Nomenclature 

Bo solar constant (W/m?) 

Oan solar declination of the day in the middle of each 
month (radians) 

dy day number of the year (dimensionless) 

Eo dm eccentricity correction factor of the earth’s orbit of 


the day in the middle of each month (dimensionless) 

Gc nm(0) monthly average hourly global irradiance on a 
horizontal surface (W/m?) 

Gg nm(O) monthly average hourly direct irradiance on a 
horizontal surface (W/m?) 

Gp nm(O) monthly average hourly diffuse irradiance on a hor- 
izontal surface (W/m?) 

Gg nm(S) monthly average hourly direct normal irradiance 
(W/m?) 


Hg am(O) monthly average daily global irradiation on a hori- 
zontal surface (Wh/m7) 

Hgo am(0) monthly average daily extraterrestrial irradiation on 
a horizontal surface (Wh/m7) 

Hg am(O) monthly average daily direct irradiation on a 
horizontal surface (Wh/m7) 

Hp am(0) monthly average daily diffuse irradiation on a 
horizontal surface (Wh/m?) 

Hg am(S) monthly average daily direct normal irradiation 


(Wh/m?) 
Hg a(S) annual direct normal irradiation (Wh/m7) 
Krm clearness index (dimensionless) 
Kpm diffuse fraction of global radiation (dimensionless) 
p latitude (°) 
w solar hour angle (radians) 
@s dm Sunrise hour angle of the day in the middle of each 


month (radians) 


are expensive and rarely available due to the cost and sophistica- 
tion of measurement devices and data processing requirements. 
There is a particular lack of data in the Sunbelt countries, which 
are more favorable for the use of CPV. 

This lack of DNI data contrasts with the availability of global 
horizontal irradiation data, which are easy to find or measure. So 
it seems to be interesting to have a procedure that obtains DNI 
from global irradiation. 

We can find in the literature several models for calculating DNI 
from different parameters [4-8]. These models are accurate but 
some of them are complex and difficult to apply, especially 
because they require input data that is not easy to obtain, for 
instance measurements taken from satellite images. In this paper, 
a simplified method obtained from previous models is presented. 
Its aim is not to improve the existing models, but to easily provide 
DNI only by using latitude and global horizontal irradiation as 
input data. The accuracy of the method seems to be acceptable for 
its use in CPV applications. 


2. DNI data sources 


Preferably, Typical Meteorological Year (TMY) data should be 
used for the analysis and design of CPV systems. These datasets 
are based on measurements taken from ground meteorological 
stations for ten years [3]. However, due to equipment costs, such 
measurements are scarce and TMY datasets are available for only 
a very few locations. An alternative solution to this problem is the 
use of a model that derives DNI from satellite data. In recent 
years, several methods for estimating direct normal irradiance 
and irradiation from satellite data have been proposed [4-8]. 
Nevertheless, the use of these procedures is often difficult as is 
obtaining the satellite data. Thus, spatial databases that directly 
provide DNI are usually utilized. 

There are a number of spatial databases that provide DNI 
values for different places and time intervals (Meteonorm [9], 
NASA-SSE [10], PVGIS [11], PVSAT [12], Satel-Light [13], SoDa 
[14], SOLEMI [15], etc.). These databases use different kinds of 
input data (satellite data, ground-based measurements) and 
different procedures in order to estimate DNI values. They are a 
very useful tool for the design and evaluation of CPV systems. 
However, some databases show weaknesses: there is a lack of 
continuity in the data because the information comes from a 
limited number of ground stations and there are no hourly DNI 
data, except for a very few sites. Although the main problem is 


uncertainty because DNI values provided by the different data- 
bases show large dispersion. The result is a fragmentation of 
services, each having its own mechanism of access and all giving 
different outputs due to different methods, input data and 
base years. 

In order to analyze this uncertainty, a comparison between 
measured DNI values taken from ground meteorological stations 
located at three Spanish places (Madrid, 40.417N -3.700W; 
Murcia, 37.983N 1.130W and Jaén, 37.766N-3.790W) and DNI 
values provided by three databases (SoDa, NASA-SSE and PVGIS) 
for these three locations has been made. 

DNI measurements for Madrid have been carried out by the 
Solar Energy Institute (Instituto de Energia Solar, IES) in 2006, 
2007 and 2010. DNI measurements for Murcia have been carried 
out by the Meteorology State Agency (Agencia Estatal de Meteor- 
ologia, AEMET) from 2005 to 2010. And DNI measurements for 
Jaén have been carried out by the University of Jaén in 2001, 2002, 
2003 and 2010. In every case, a pyrheliometer with wavelength 
range from 0.2 to 4 um has been used for the DNI measurements. 

SoDa website disseminates the HelioClim (HC) databases 
which contain radiation values at ground level obtained from 
the processing of the images taken by Meteosat satellites. In this 
case, the HC databases are processed by MINES ParisTech. On the 
other hand, NASA-SSE release 6 combines results from GEWEX/ 
SRB 3 and ISCCP projects with NCAR reanalysis products (1983- 
2005). The database directly provides monthly and annual DNI 
values around the world. Finally, PVGIS includes solar radiation 
developed by a combination of solar radiation models and 
interpolated ground observations, representing the period 
1981-1990. It does not directly provide DNI but this can be 
estimated from global and diffuse radiation components. 

Table 1 compares measured values taken from the meteor- 
ological stations to values provided by the databases. Annual DNI 
value measured at Madrid is 1882 kW h/m? while the different 
databases give values between 1814 and 1964 kW h/m?. This 
implies an error between —3.6% and 4.4%. Annual DNI value 


Table 1 
Annual DNI values for three Spanish locations (kW h/m7). 


Location Measured SoDa NASA-SSE PVGIS 
MADRID 1882 1947 1964 1814 
MURCIA 1961 1738 1905 1870 
JAEN 1984 2208 2132 1745 
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measured at Murcia is 1961 kW h/m? while the different data- 
bases give values between 1738 and 1905 kW h/m7. This implies 
an error between —11.4% and —2.9%. And annual DNI value 
measured at Jaén is 1984 kW h/m? while the different databases 
give values between 1745 and 2208 kW h/m?. This implies an 
error between —12.0% and 11.3%. Fig. 1 shows the annual DNI 
error for the three databases. A clear trend cannot be observed, 
although PVGIS underestimates DNI for these three locations and 
the three databases underestimate Murcias DNI value. As a 
conclusion, it can be assumed that it is difficult to get reliable 
DNI data because each database gives different values and thus, 
there are sources of uncertainty. 

Other authors have reported the same conclusions about uncer- 
tainty [16,17]. Sari et al. [16] cross-compared maps of annual DNI 
from five databases offering solar resource and climate data for 
Europe: Meteonorm [18], NASA-SSE [19], PVGIS [20], Satel-Light [21] 
and SOLEMI [22]. From these databases, the following maps were 
calculated: (1) overall average estimated by a simple averaging, 
(2) five maps of differences to the average, and (3) standard devia- 
tion. While maps of differences show the deviation of the particular 
dataset from the overall average, the standard deviation indicates 
users uncertainty given by the magnitude of differences between the 
combined data sources. Results show that 90% of the analyzed area 
presents an uncertainty in annual DNI, expressed by relative stan- 
dard deviation, of up to 17%. This high uncertainty indicates that it is 
necessary to improve the databases and the radiation models. Some 
efforts have been made in this field. For instance, the MESoR project 
[17], funded by the European Commission, reduces the associated 
uncertainty by setting up standard benchmarking rules and mea- 
sures for comparing the databases, user guidance to the application 
of resource data and unifying access to various databases. 

Below, a very simple procedure for the calculation of DNI 
based on previous existing models is presented. It can be easily 
implemented in a spreadsheet or in computer applications in 
renewable energy. The aim of the proposed procedure is not to 
provide high accuracy. It tries to meet the basic requirements for 
the analysis of CPV systems by providing DNI errors similar to 
those of the different meteorological databases. Inputs of the 
procedure are the global horizontal irradiation values which, 
unlike DNI values, are often measured in the meteorological 
stations and are easy to find. The procedure does not improve 
the existing methods but is an easy-to-use choice useful for 
teaching and providing acceptable accuracy. 


3. Description of the proposed method 


In this section, a procedure that calculates DNI values from 
global horizontal irradiation values is presented. 


15% 
10% 


5% 
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Fig. 1. Annual DNI error for the selected locations and databases. 


In what follows, according to [23,24], letter G will be used to 
indicate irradiance or the power incident on a surface per unit 
area of surface (units: W/m?) and letter H will be used to indicate 
irradiation or the energy incident on a surface over a specified 
period of time per unit area of surface (units: W h/m?). 

At the high concentration rates involved in High Concentration 
Photovoltaic (HCPV) systems, time resolution as well as uncer- 
tainty in radiation data play a complex role [25]. However, some 
authors have proven that the complexity of a solar radiation 
model needed for yearly energy calculations is very low [26]. 
Twelve values of monthly mean solar irradiation are enough to 
estimate yearly energy with small errors. Also, time resolutions 
shorter than hourly samples do not significantly improve the 
energy estimation results. 

So, the start-up data required for the calculations are the site 
latitude, ~, and the twelve values of monthly average daily global 
irradiation on a horizontal surface, H¢gm(0). Results of the method 
are the monthly average direct normal irradiance at one hour 
intervals for the twelve months of the year, Gg nm(S), the twelve 
values of monthly average daily direct normal irradiation, Hg gm(S), 
and the value of the annual direct normal irradiation, Hg a(S). 
The proposed method performs the sequence shown in Fig. 2 in 
order to calculate the DNI values. 

The method is based on easy-to-use equations proposed by 
several authors [27] which provide a good approximation for a wide 
range of latitudes. Other more accurate equations and correlations 
[28-30] have not been taken into account because of their complex- 
ity or because they are only valid for a specific area. 

The proposed method is based on the following procedure: 


I. Solar declination of the day in the middle of each month, in 
radians, ĉan [31]: 


Oan = er sin (27 x 


dn +284 
= ) (1) 


365 


where d, is the day number of the year, ranging from 1 on 
January 1 to 365 on December 31. The values corresponding 
to January 15, February 15, March 15, etc. will be used. 

II. Eccentricity correction factor of the earth’s orbit of the day in 
the middle of each month, dimensionless, £o am [32]: 


Eo dm = 1+0.033 x cos (27 x - ) (2) 


III. Sunrise hour angle of the day in the middle of each month, in 
radians, Ws; am [33]: 


sam = —cos~'(—tan dgn x tan ~) (3) 


1. Calculation of daily extra terrestrial irradiation ona horizontal surface. 
2. Calculation of daily diffuse irradiation on a horizontal surface. 

3. Diffuse irradiance calculation from irradiation. 

4. Global irradiance calculation from irradiation. 

5. Calculation of hourly direct normal irradiance. 

6. Calculation of daily direct normal irradiation. 

7. Calculation of annual direct normal irradiation. 


Fig. 2. Flow diagram of the proposed method. 
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IV. Monthly average daily extraterrestrial irradiation on a hor- 
izontal surface of each month, in W h/m?, Hgo am(0) [34]: 


24 
Hgo dm(0) = — Bo x €0 am x (COS : COSÒan) 

X(Dsam X COSMsdm—SIN@ sm) (4) 
being Byp=1367 W/m? the solar constant [35]. It is assumed 
that the extraterrestrial irradiation of the day in the middle 
of each month almost equals the monthly average daily 


extraterrestrial irradiation of this month. 
V. Monthly average clearness index, dimensionless, Krm [36]: 


Krm = Hgam(0)/Hgoam(0) (5) 


VI. Monthly average diffuse fraction of global radiation, dimen- 
sionless, Kpm [37]: 


Kpm = 1-1.13Ktm (6) 


VII. Monthly average daily diffuse irradiation on a horizontal 
surface of each month, in Wh/m?, Hp am(0): 


Hp am(0) = Hgam(0) x Kpm (7) 


VII. Monthly average daily direct irradiation on a horizontal 
surface of each month, in W h/m?2, Hg am(0): 


Hg dm(0) = Hg dm(0)—Hp am(0) (8) 


IX. Monthly average hourly diffuse irradiance on a horizontal 
surface, in W/m?, Gp nm(0) [38]: 


T COS@—COS(. 
Gp tm() = Hp dm(0) 53 * ( A ) 


Ms dm X COS@sdm— SİNO; dm 


(9) 


where œ is the solar hour angle in radians (—z for 0:00 
hours and +2 for 24:00 hours). Irradiance values are 
calculated at one hour intervals (A@=27/24 rad). 

X. Monthly average hourly global irradiance on a horizontal 
surface, in W/m?, Ge nm(0) [39]: 


T 
Gg hm(0) = Hg am(0) 576 +b 


7 cosc) ( COSW—COSWsdm ) (10) 
Os dm X COSMsdm—SIN®s dm 
where a and b parameters are calculated by the formulae: 
a= 0.4090—0.5016 sin(@, am + 1.047) (11) 
b = 0.6609 + 0.4767sin(@s am + 1.047) (12) 


XI. Monthly average hourly direct irradiance on a horizontal 
surface, in W/m’, Gg nm(0): 


Gg hm(0) = Gg hm(0)—Gp hm(0) (13) 


XII. Monthly average hourly direct normal irradiance, in W/m?, 
Gg hm(S) [40]: 


Gg hm(0) 
sindgn X SIN@+COSdgn X COSP x COSM 


Germ (S) = (14) 


XIII. Monthly average daily direct normal irradiation, in Wh/m?, 
Hg dm(S): 


24 
Ag am(S) © 5 GB nmj(S) x At (15) 
j=l 


where Gg nm j(S) is the monthly average hourly direct normal 
irradiance at “j” hour and At=1 h. 


XIV. Annual direct normal irradiation, in W h/m?, Hg a(S): 


12 
Hs a(S) ~= X` dm; x Hg am i(S) (16) 


i=1 


where dm; is the number of days for the “i” month and Hg am £S) 
is the monthly average daily direct normal irradiation for the 


um 


i” month. 


4. Results and discussion 


In this section, the proposed method is applied to three 
Spanish locations (Madrid, Murcia and Jaén) in order to estimate 
both the monthly average daily direct normal irradiation of each 
month, Hg am(S), and the annual direct normal irradiation, Hg a(S). 
These locations correspond to the locations analyzed in Section 2. 
They have been selected because ground-based measurements of 
annual DNI are available and thus, a comparison between mea- 
sured and calculated values can be made. In what follows, 
irradiation values are expressed in kW h/m’. 

Table 2 shows the input data used for the calculations, i.e., the 
site latitude, @, and the twelve values of monthly average daily 
global irradiation on a horizontal surface, Hc am(0). The global 
irradiation values have been obtained from PVGIS. 

With these input data, the procedure described in Section 3 can 
be easily applied. Results are shown in Table 3, including the twelve 
values of monthly average daily direct normal irradiation, Hg am(S), 
and the value of the annual direct normal irradiation, Hg a(S). 

In order to compare measured and calculated DNI values, the 
annual DNI error has been computed for each location: 


Hp alS)catcutated HB a(S)measured 400 (17) 
Hg a(S) measured 


Table 4 shows the results. The proposed method overestimates 
DNI for the three locations. Madrid’s error is 3.9% while the three 
meteorological databases selected in Section 2 (SoDa, NASA-SSE 
and PVGIS) gave errors between — 3.6% and 4.4%. Murcias error is 
1.2% while the databases errors were between —11.4% and 
—2.9%. And Jaéns error is 0.2% while the databases errors were 
between — 12.0% and 11.3%. As a conclusion, annual DNI values 
produced by the proposed method are comparable to those of the 
existing meteorological databases. Thus, an alternative easy-to- 
use method can be used to compute DNI. 


error(%) = 


Table 2 
Input data required by the proposed method for the three selected locations: 
latitude and monthly average daily global irradiation on a horizontal surface. 


Location: Madrid Murcia Jaén 
latitude, °: 40.417 37.983 37.766 
Month He am(0), KW h/m? 

J 1.98 2.44 2.46 
F 2.69 3.21 3.16 
M 4.44 4.54 4.65 
A 5.09 5.40 5.29 
M 6.50 6.48 6.59 
J 7.23 7.00 7.19 
J 7.33 7.05 7.07 
A 6.44 6.15 6.22 
S 4.98 5.00 4.94 
(0) 3.36 3.74 3.74 
N 2.14 2.47 2.49 
D 1.61 2.16 2.09 


PJ. Pérez-Higueras et al. / Renewable and Sustainable Energy Reviews 16 (2012) 5529-5534 5533 


5. Application 


The simplified procedure described in this paper can easily be 
used to obtain DNI data worldwide. The only requirement is to 
know global horizontal irradiation values for the area of interest. 
As these values are easy to find or measure, the proposed method 


Table 3 
Output data obtained by the proposed method: monthly average daily DNI of each 
month and annual DNI. 


Location: Madrid Murcia Jaén 
Month Hp am(S), kW h/m? 

J 3.12 4.02 4.02 
F 3.44 4.29 4.10 
M 5.72 5.44 5.69 
A 5.04 5.39 5.14 
M 6.75 6.55 6.78 
J 7.86 7.19 7.60 
J 8.37 7.61 7.65 
A 7.55 6.58 6.72 
S 6.20 5.79 5.62 
(0) 4.56 5.03 4.98 
N 3.21 3.64 3.65 
D 2.44 3:73 3.42 

Hg a(S), KW h/m? 
Annual 1955 1985 1988 
Table 4 


Comparison between measured annual DNI and annual DNI calculated by the 
proposed method. 


Location Hpg a(S), KW h/m? 

Measured Calculated Error (%) 
Madrid 1882 1955 3.9 
Murcia 1961 1985 1.2 
Jaén 1984 1988 0.2 


becomes a very useful tool in renewable energy applications, 
particularly in CPV applications. 

As an example of such applications, a complete map of the 
Spanish annual DNI is presented in Fig. 3. It has been obtained 
from global irradiation data of the 47 peninsular provincial 
capitals of Spain. For each location, annual DNI has been calcu- 
lated by using a Microsoft Excel spreadsheet with the procedure 
described. This spreadsheet is available at http://www.ujaen.es/ 
investiga/solar/dni.xlsx. The DNI data obtained has been graphi- 
cally represented with the help of the Surfer program by Golden 
Software [41]. 

This example proves that the proposed method can be used to 
obtain a mass of DNI data for different regions worldwide with 
acceptable accuracy. 


6. Conclusions 


DNI data provided by the existing meteorological databases is 
obtained from satellite images and ground-based measurements 
by using different methods, input data and base years. While 
these databases are a useful tool in renewable energy applica- 
tions, they present uncertainty. 

A simplified procedure for estimating DNI from global irradia- 
tion based on previous existing models has been developed. It 
provides accuracy similar to that of the existing databases. As 
global irradiation data are easy to find or measure, the method 
becomes a good choice for the design and analysis of CPV and 
other renewable energy systems. It can be easily implemented in 
computer programs and is useful for teaching. 

The accuracy of the method has been shown by comparing 
calculated and measured DNI for three Spanish locations. Also, the 
estimation error has been compared with the errors provided by 
three existing databases. Results indicate that the procedure can 
be used as an alternative to other sources of DNI data. 

In order to prove the ability of the proposed method to easily 
obtain a mass of DNI data for different areas worldwide, a map of 
the Spanish annual DNI has been developed and presented. 


Fig. 3. Map of the Spanish annual DNI provided by the proposed method. 
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